Trapped Atomic Fermi Gases 
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Abstract 

A many-body system of fermion atoms with a model interaction characterized by the scattering 
length a is considered. We treat both a and the density as parameters assuming that the system can 
be created artificially in a trap. If a is negative the system becomes strongly correlated at densities 
p ~ | a | 3 , provided the scattering length is the dominant parameter of the problem. It means that 
we consider \a\ to be much bigger than the radius of the interaction or any other relevant parameter 
of the system. The density p c \ at which the compressibility vanishes is defined by p c \ ~ | ct | 3 . Thus, 
a system composed of fermion atoms with the scattering length a — > — oo is completely unstable at 
low densities, inevitably collapsing until the repulsive core stops the density growth. As a result, any 
Fermi system possesses the equilibrium density and energy if the bare particle-particle interaction 
is sufficiently strong to make a negative and to be the dominant parameter. This behavior can be 
realized in a trap. Our results show that a low density neutron matter can have the equilibrium 
density. 
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Although a theory of Fermi gases is not yet well developed, there are nevertheless new challenging 
experimental possibilities to explore trapped Fermi gases 0. These are expected to stimulate a proper 
theoretical description of these Fermi quantum systems. Currently this description of quantum Fermi 
gases at their different states, which can be equilibrium, quasi-equilibrium or far from equilibrium, is 
usually based on tedious numerical calculations, particularly when the interaction between particles, or 
O . atoms, of a gas cannot be considered as weak. This implies that the dimensionless effective coupling 
constant of the interaction ppa > 1. Here pp is the Fermi momentum and is related to the system's 
density by p = pp/37i 2 , while a is the scattering length. There are a few cases, when it is possible to treat 
the properties analytically. The Random Phase Approximation (RPA) is applicable for a high density 
electron gas and the low density approximation deals with dilute gases ||. In both cases the kinetic 
energy Tj, is assumed to be much bigger than the interaction energy E int of the system. This permits the 
application of some kind of a perturbation theory. In the case of an electron liquid it turns out that the 
analytical RPA-like description is also possible not only at very high but also at medium densities when 
^ .Xjfe ~ Em @, |1- Similar extension of the range of validity is impossible in the case of fermion systems 
at low densities p (there the gas approximation is not applicable if ppa > 1, or Tj. ~ E int ). If the pair 
interaction is attractive and sufficiently strong, the system can have a quasi-equilibrium or equilibrium 
states in which T k ~ E int . On the other hand, these states are separated by special regions at which the 
incompressibility K(p) < and the system is completely unstable. An experimental study performed on 
such Fermi-systems would be of great importance presenting new information on the behavior of dilute 
gases and on the gas-liquid phase transition. The observation of the existence of such regions and points 
at which K(p) = can present a challenging problem for a theory designed to describe these peculiarities. 
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In this Letter we address the above mentioned problem and consider an infinitely extended system 
composed of Fermi particles, or atoms, interacting by an artificially constructed potential with the desirable 
scattering length. We demonstrate that the consideration can be accomplished analytically provided that 
the pair interaction between fermions is characterized only by the scattering length. That is \a\ is much 
larger than the other relevant parameters of the interaction, or a — ► — oo. In this case one can say that the 
scattering length is the dominant parameter of the problem under consideration. As it will be demonstrated 
below, in such a case the system is located in the vicinity of the gas-liquid phase transition, transforming 
it into a strongly correlated one. Therefore, the problem of calculating its properties has to be treated for 
the most part qualitatively. Such an investigation is of great importance since it can be applied to fermion 
systems interacting via potentials with not only infinite, but also sufficiently large a. For instance, the 
scattering length a of neutron-neutron interaction is about —20 fm, thus greatly surpassing the radius of 
the interaction tq. Then this investigation can be viewed as the first step to study trapped Fermi gases, 
which are systems composed of Fermi atoms interacting by an artificially constructed potential with almost 
any desirable scattering length, similarly to that done for the trapped Bose gases, see e.g. 0. 

Let us start by considering the general properties of a Fermi system with some attractive two-particle 
bare interaction gV(r) of range r and g defining the strength of the interaction. We assume that gV(r) 
is sufficiently weak to create a two-particle bound state and that the scattering length a corresponding 
to this potential is the dominant parameter, being negative and finite, while r /|a| <C 1. Then in the 
Hartree-Fock approximation the ground state energy Ehf{p) is given by the following expression 

E HF (p) = ^p + t HF p 2 , (1) 

where M is the particle mass, and the parameter tjjp , being negative, is entirely determined by the 
potential gV(r). For instance, in the case of a short range 5-type interaction one has tuF = ~g/^- Eq. 
(1) shows that at low densities E HF > due to the kinetic energy term, but at sufficiently high densities 
p — > oo the potential energy tuFP 2 becomes dominating, leading to the collapse of the system, with 
Ehf — > — oo. Keeping in mind that the Hartree-Fock approximation gives the upper limit to the binding 
energy Ehf > E, one can conclude that the system does not have, in this case, an equilibrium density 
p e and energy E e since E e — > — oo when p — > oo ||. Note, that for a given and finite total number of 
particles N, the HF energy is not going to infinity and the system collapses into a small volume with the 
radius r , with the density p ~ N/r^. At the densities p — > we deal with a dilute gas and the energy 
E(p) — > 0, remaining positive at these densities J7|. Therefore, it must have at least one maximum at the 
density p m before it becomes negative, on the way to E — > — oo. If the potential gV(r) includes a kind of 
a "repulsive core" at sufficiently short distances, the system has an equilibrium density and energy, p e and 
E e , respectively, determined by the repulsive core strength and its radius r c ~ r . The general features of 
the function E(p) can be qualitatively represented by a simple expression 

E(p) = -^p + t p 2 + t 3 p 3 . (2) 

The first term of Eq. (2) is the kinetic energy Tf., while the second and the third terms are related to the 
interaction energy E int determined by the potential gV(r). The second term which is proportional to to 
gives a qualitative description in the gas limit. The third term provides the behavior of E(p) at higher 
densities, including that of the equilibrium density, so that t < 0, and < t 3 . 

Now we apply Eq. (2) to demonstrate the most important features of the system under consideration: 
a) when p — > the third term on r.h.s. in Eq. (2) can be omitted. The kinetic energy is relatively 
very big, T k 3> E int , and t ~ a, with a < being the scattering length. In that case we have a dilute 
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Fermi gas with positive pressure P and incompressibility K, the latter being determined by the equation, 
see e.g. [g, 

2 dE 2 {p) 

K (P) = P TJ> ■ ( 3 ) 



b) on the way to higher densities, which can be achieved by applying an external pressure, the system 
reaches the density p cl < p m at which the incompressibility is equal to zero, K(p cl ) = 0. Remembering 
that at the maximum the second derivative is negative, one can conclude, as it is seen from Eq. (3), that 
K(p m ) < 0. In the range p c2 > p > p c \ the incompressibility is negative, K < 0, and as a result the 
system becomes totally unstable. In fact, in this density range all calculations of the ground state energy 
are meaningless since such a system cannot exist and thus be observed experimentally ||; 

c) at some point p = p C 2 > Pm the contribution due to the repulsive core becomes sufficiently strong to 
prevent the further collapse of the system. The incompressibility attains K = at p C 2 < p e , being positive 
at the higher densities. Finally, the system becomes stable at p > p c2 , reaching equilibrium density at 
p e with equilibrium energy equal to E e . Generally speaking, the system is quasi-stable at the densities 
< p < pd and pci < p < p e because, being affected by a sufficiently strong external field, it occupies 
its real stable state with the energy E e and density p e 0. For the sake of simplicity, we call these states 
stable states. It is obvious that K(p e ) > being proportional to the second derivative at the minimum [see 
Eq. (3)]. It should be kept in mind that in this density domain, p > p c2 , the function E(p) is determined 
by the repulsive part of the potential which makes t% > 0. As mentioned above, without this component 
of gV{r) the system's energy would decrease, E HF (p) — > — oo, with the density growth, p — > oo, thus 
inevitably collapsing. Indeed, in this case, the potential being pure attractive has no structure to ensure 
any stable states at the densities p > p c \. As a result, one can write down a dimensionless expression for 
the ground state energy as a function of the only variable z = pfo> [0, 

aE(z)=z 5 (l + /3(z)) : (4) 

with a = l0n 2 Ma 5 . Note, that the function (3{z) dependents only the variable z. In the low density limit, 
\app\ <C 1 and when the interaction has the radius r , Eq. (4) reads M, 



aE(z) = z 
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Here the function 7(7/, z) is of the order of one, 7(^,2;) ~ 1. It is seen from Eq. (5) that as soon as 
the scattering length becomes large enough, \a\ 3> ro, one can omit the contribution coming from the 
function 7 and neglect all the term proportional to (r /a) 3 . Then Eq. (5) reduces to Eq. (4). Thus, in 
the case \a\ — » 00 we can use Eq. (4) to determine the ground state energy E. Eq. (4) is valid up to the 
density p c \ which is a singular point of the function ft(z), since beyond this point K < 0, and the system 
is completely unstable. On the other hand, there is no physical reason to have another irregular point 
in the region < p < p c \. Using Eq. (3) for the incompressibility and Eq. (4) for the energy, one can 
calculate the position of the point z c \ where K = 0. Denoting the corresponding z as z c \ = Co , where Co 
is a dimensionless number, one is led to the conclusion that p c \ ~ |a|~ 3 provided a is sufficiently large to 
be the only dominating parameter. As it is seen from Eq. (5), we can expect corrections of the order of 
(ro/a) 9 to this universal behavior. Thus, the system has only one stable region at small densities p < p c \ 
which decreases and even vanishes as soon as a — > —00, and the function (3 is in fact determined only in 
the region \z\ < \z&\ |IIJ. 



Consider the behavior of the system when the density approaches p — > p c \ ~ \a\ from the low density 
side. Normally, points p c \ and p c2 , overlooked in calculations because of the lack of self consistency 
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which relates the linear response function of the system to its incompressibility K, 



X(q -> 0,iu -> 0) 



'd 2 E s 



These points can give important contributions to the ground state energy, 
energy of the system in the following form (see e.g. flTTf), 



(6) 



To see this we express the 



E(p) = E HF (p) - - J \x(q, iu, g) - Xo(q, «*>)] v(q) 



dq du) dg 



(7) 



where xiq^^^d) is the linear response function on the imaginary axis, xoiq,^) is the linear response 
function of noninteracting particles, and v (q) is the Fourier image of gV(r). The integration over uj goes 
from — oo to +oo, while the integration over the coupling constant g runs from zero to the real value of the 
coupling constant, i.e. to g = 1. At the point p = p cl the linear response function has a pole at the origin 
of coordinates q = 0, to = due to Eq. (6). At the densities p > p c \ the function x(q, iv,9) has poles at 
finite values of the momentum q and frequencies iu. This prevents the integration over iu, making the 
integral in Eq. (7) divergent. Thus, we conclude that it is the contribution of these poles that reflects the 
system's instability in the density range p c \ < p < p c2 . Note, that violations of Eq. (6) lead to serious 
errors in the calculation of the ground state energy. Equation (7) can be rewritten, explicitly accounting 
for the effective interparticle interaction R(q,iu, g), in the following form 



E(p) 



Ehf(p) ~ - 



and x is given by the following equation 

x(q,u,g) 



Xo(q,iu>)R(q,iu,g) 
1 - R(g,iu,g)xo(q,iw) 

m 

xo(g,w) 



v(q) 



dq duj dg 
2(2tt) 4 : 



(8) 



R{q,^,g)xo(q,u;) 



It is seen from Eqs. (6) and (9) that the denominator (1 — Rxo) vanishes at p 
obtain 

1 1 

R(q -> 0, u = 0, g = 1) oc 

Xo(g,0) p F M + q 2 /p z F 

Here pp ~ p c \. If the scattering length tends to infinity driving p cl — > 0, it follows from Eq. (10) that 



(9) 

p c \. As a result, we 
(10) 



R(q 



0,0,1) a 



(11) 



and we can conclude that the effective interaction R behaves as a gravitational-like field, leading the 
system to collapse until the repulsive core stops the further squeezing of the system. Note, it is impossible 
to represent the denominator as a power series in Rxo through approximating the expansion by a finite 
number of terms. This result is quite obvious since p c \ is a branch point in the function E(p), which makes 
it impossible to expand that function in the vicinity of this point. Therefore, one should try to satisfy Eq. 
(6) in order to get proper results for the ground state calculations in the vicinity of the instability points 

§B0 



We would like to demonstrate here the essential qualitative features of atoms in a trap, which can be 
observed experimentally at the densities near the critical density, p < p c \. It follows from Eq. (9), at least 
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for small q and p — > p c \, that the susceptibility x(q 0, 0, 1), being proportional to the incompressibility 
(see Eq. (6)), is divergent. Therefore, the denominator in Eq. (9) can be expressed in the form 

1 - R(q, 0, l) X o(q, 0) ~ a {p/p el ) + a ig 2 /p 2 F . (12) 

Here «o(p/p c i) is a positive function such that a (l) = 0, and «i is a positive constant. As a result, the 
linear response function takes the form, 

X (g-> 0,0,1)1^ oc^. (13) 

Since the linear response function is the density-density correlation function, we obtain from Eq. (13) 
that the radius of the correlation tends to infinity as it must in the vicinity of the phase transition point 
||. For instance, the density fluctuations 5p(q) induced by an external potential, Sp(q) = x(<l, 0)v ext (q) oc 
v ext (q)/q 2 . Thus, an external field v ext (q) even it is of a short range cannot be screened and produces the 
fluctuations of the density of extremely long range. It seems that such a behavior can be realized in a trap 
through the observation of a strong light scattering by the density fluctuations of multi-fermion system 
when p —>■ p c \. On the other hand, the same picture can be observed if the density is fixed at some po value, 
while the critical density is driven to this density by adjusting the scattering length by an external field. 
If the bare potential were purely attractive, the interval of the densities [0,p c i] within which the system is 
stable, would vanish with the growth of \a\ since p c \ ~ M~ 3 . As a result, in the limit a = — oo the density 
Pd — > 0, and the considered system would be completely unstable at any density. As soon as the scattering 
length deviates from its infinite value, that is +oo > a > — oo the system comes back to its stable state 
at list in the range of the density values p < p c \ ~ |a|~ 3 . Note, as it follows from our consideration, 
that any Fermi system possesses an equilibrium density and energy if the bare particle-particle interaction 
contains a repulsive core and its attractive part is strong enough, so that a — > — oo. Indeed, at sufficiently 
small densities the energy is negative. The system collapses (since at these densities the incompressibility 
K < 0) until the core stops the density growth. Therefore, the minimal value of the ground state energy 
must be negative when the repulsive core switches on to prevent the system from further collapse. 

Note that trapped alkali Fermi gases are spin polarized and cannot interact in s-wave channel. Consider 
a two component system realized by mixtures of two Fermi alkali gases, say 6 Li and 40 K. We assume that 
the corresponding number densities p\ and p 2 of these gases coincide, p\ = p 2 = p. In that case, we 
come back to our case: at each level there are two spin-polarized atoms interacting in s-wave and with 
"effective" mass M* = M\M%/(Mi + M 2 ). Here M\ and M 2 are the particle masses of the alkali gases 
under consideration. As a result upon using Eq. (4), we obtain p cl ~ |a c |~ 3 , where a c is the scattering 
length related to the interaction between the different alkali atoms. It is possible to consider a large variety 
of two component Fermi systems with p% ^ p 2 , or even two component systems composed of Fermi and 
Bose gases which can retain the significant features of the considered one component Fermi system. A 
detailed investigation of such systems is in progress and will be published elsewhere. It is worth remarking, 
that superfluid correlation cannot stop the system from squeezing, since their contribution to the ground 
state energy is negative. After all, the superfluid correlation can be considered as additional degrees of 
freedom which can therefore only decrease the energy. Recent considerations show that the compressibility 
related to the superfluid correlation is positive (see, e.g. jn| and references therein). Therefore we do not 
expect that the contribution coming from the superfluid correlation can lead to a significant change in our 
estimate of the critical density, p c \ ~ |a|~ 3 . While deriving the exact relationship, one has to include a 
proper treatment of the pairing. 

A liquid similar to the model considered above exists in Nature, viz. liquid 3 He. If a 3 He2 dimer exists, 
obviously its bounding energy does not exceed the bounding energy of a 4 He2 dimer which is 10 -4 meV WA 



The ground state energy of liquid helium is about 2*10 meV per atom. Because of this huge difference 
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in binding energies, it is evident that the contribution coming from the binding energy of the dimer to the 
ground state energy of the liquid is insignificant. In fact, numerical calculations show that the pair potential 
is rather weak to produce the dimer 3 He2 Hl5| . Thus, one can reliably consider an infinite homogeneous 



system of Helium atoms as consisting of particles interacting via pair potentials, characterized by a very 
large but finite scattering length \a\ ^> r^. The following additional remark is appropriate. It seems 
quite probable that the neutron-neutron scattering length (a ~ —20 fm) is sufficiently large to permit the 



neutron matter to have an equilibrium energy and density [12]. 

In summary, a system of fermions interacting by an artificial potential has been considered. The 
qualitative consideration presented above gives strong evidence that the system becomes unstable at 
densities p > p c \. In the vicinity of p c \ ~ H~ 3 , there are long-ranged fluctuations of the density which 
can be realized in a trap through the strong scattering of light. Our results suggest that the equation of 
state of a low density neutron matter has peculiarities. 
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